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ABSTRACT 

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and 
the Galaxy and Mass Assembly project, we show that, for galaxy masses above 
~ 10 8 M 0 , 51% of the stellar mass-density in the local Universe is in early-type galax¬ 
ies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density 
is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift 
benchmark, we have used a calorimetric technique to quantify the importance of the 
morphological transformation of galaxies over the history of the Universe. The extra- 
galactic background radiation contains all the energy generated by nuclear fusion in 
stars since the Big Bang. By resolving this background radiation into individual galax¬ 
ies using the deepest far-infrared survey with the Herschel Space Observatory and a 
deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using 
measurements of the Sersic index of these galaxies derived from the HST images, we 
estimate that ~83% of the stellar mass-density formed over the history of the Universe 
occurred in LTGs. The difference between this and the fraction of the stellar mass- 
density that is in LTGs today implies there must have been a major transformation 
of LTGs into ETGs after the formation of most of the stars. 

Key words: galaxies: bulges, evolution, star formation 


1 INTRODUCTION 

Over the last decade much of the research on galaxy evo¬ 
lution has started from the paradigm that there are two 
types of galaxy: star-forming galaxies and passive or quies¬ 
cent galaxies in which stars are no longer forming. In this 
paradigm star-forming galaxies lie on the Galaxy Main Se¬ 
quence (GMS) until some catastrophic process (or at least, a 
process brief compared to the age of the Universe) quenches 
the star formation, causing the galaxy to move away from 
the GMS and become a ‘red and dead’ galaxy (e.g. Peng et 
al. 2010; Speagle et al. 2014). 

The rapid evolution of the stellar mass function of pas¬ 
sive galaxies since z ~ 1 — 2, relative to the much slower 
evolution of the stellar mass function of star-forming galax¬ 
ies (Faber et al. 2007; Muzzin et al. 2013; Ilbert et al. 2013), 
implies that the formation of passive galaxies occurred after 
the epoch at which the star-formation rate in the Universe 
was at its peak (Hopkins and Beacom 2006). This tempo¬ 
ral sequence is circumstantial evidence that passive galax¬ 
ies were formed from the star-forming population. A more 
quantitative argument is that Peng et al. (2010) have shown 
that the shapes of the two stellar mass functions are exactly 
what one would predict from a simple model in which passive 
galaxies are formed by the quenching of the star formation 
in the galaxies in the star-forming population. 

In an accompanying paper (Eales et al. 2015), we show 
that when this paradigm is considered in the light of the 
galaxies detected by the Herschel Space Observatory (Pil- 
bratt et al. 2010) it looks less convincing. Approximately 
25% of the galaxies detected at low redshift by the largest 
Herschel extragalactic survey, the Herschel ATLAS (Eales 
et al. 2010), are galaxies that would have been classified as 
quiescent or passive using optical criteria, but still have large 
reservoirs of interstellar material and are still forming stars. 
They are red but not dead. The blurring of this dichotomy 
when viewed from a far-infrared perspective is, we argue, 
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evidence that some more gradual physical processes (rather 
than catastrophic quenching) are necessary to explain the 
properties of the galaxy population (Eales et al. 2015). At 
the very least, our results show that there is a major prac¬ 
tical problem in measuring star-formation rates in galaxies 
and in separating galaxies into the two classes, which may 
explain the huge diversity in measurements of the shape and 
evolution of the GMS found by Speagle et al. (2014) in their 
meta-analysis of 25 separate studies of the GMS. 

In this paper we consider a different dichotomy: that be¬ 
tween early-type and late-type galaxies (henceforth ETGs 
and LTGs). The terms ‘early type’ and ‘late type’ were 
coined by Hubble (Hubble 1926, 1927), although he did not 
mean the terminology to indicate that ETGs evolve into 
LTGs (Baldry 2008). ETGs are galaxies with a prominent 
spheroid or bulge, in particular classes E and SO in the Hub¬ 
ble morphological classification, whereas the structures of 
LTGs are dominated by a spiral pattern. We discuss below 
how to turn this rather vague definition into one that is 
more quantitative. Although there is an obvious dichotomy 
between elliptical galaxies and late-type spirals when looking 
at optical images, the dichotomy between ETGs and LTGs 
is much less obvious when galaxies are studied more closely. 
Integral-field kinematic studies have shown that ETGs can 
be split into slow and fast rotators, with the kinematic prop¬ 
erties of the fast rotators being similar to those of LTGs and 
the slow rotators being much rarer than the fast rotators 
(Emsellem et al. 2011). Careful quantitative analysis of the 
images of ETGs also often reveals that the galaxy has a disk 
even when this is not apparent from simply looking at the 
image (Krajnovic et al. 2013). Observations with the Her¬ 
schel Space Observatory have revealed that ~50% of ETGs 
contain dust, and that there is gradual change of the dust 
properties of galaxies as one moves along the Hubble Se¬ 
quence from LTGs to ETGs, rather than a sudden jump at 
the transition from late types to early types (Smith et al. 
2012a). 

Although Hubble did not draw any conclusions about 
the evolutionary links between ETGs and LTGs, there are 
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many theoretical arguments for why galaxy evolution is ex¬ 
pected to include the morphological transformation of galax¬ 
ies. Almost four decades ago, Toomre (1977) suggested that 
the merger of two LTGs might scramble the stars and ve¬ 
locity fields, leading to the formation of an ETG supported 
against gravity by the random velocities of the stars rather 
than by the rotation of a galactic disk. On the other hand, 
if one starts with a spheroid, it seems likely that the accre¬ 
tion of gas and angular momentum from the intergalactic 
medium is likely to lead to the growth of a disk (Combes 
et al. 2014 and references therein). A recent idea is that 
a spheroid is formed by the rapid motion of star-forming 
clumps towards the centre of a disk (Noguchi 1999; Bour- 
naud et al. 2007; Genzel et al. 2011, 2014; Dekel and Burkert 
2014), transforming a LTG into an ETG, with the growth 
of the spheroid eventually stopping the formation of stars in 
the disk once the spheroid has enough mass to stabilize the 
disk against gravitational collapse (Martig et al. 2009). 

There is substantial overlap between ETGs and the 
class of passive galaxies and between LTGs and the class 
of star-forming galaxies. The overlap between the classes in 
the low-redshift Universe is well known (Kennicutt 1998), 
but over the last decade a large number of observational 
programmes with the Hubble Space Telescope (HST) have 
shown that this is also true at high redshift: high-redshift 
passive galaxies tend to have spheroidal structures and the 
structures of high-redshift star-forming galaxies are gener¬ 
ally dominated by disks (Bell et al. 2004a; Wuyts et al. 2011; 
Bell et al. 2012; Bruce et al. 2012; Buitrago et al. 2013; 
Szomoru et al. 2013; Tasca et al. 2014). These studies show 
that the evolutionary process that produces passive galaxies 
must be linked to the evolutionary process that generates 
spheroids (Bell et al. 2012). 

These HST studies have two limitations. The first lim¬ 
itation is a consequence of the extragalactic background ra¬ 
diation and of the nature of the sources which constittue 
it. The extragalactic background radiation (Dole et al. 
2006; Dwek and Krennrich 2013) is dominated by radia¬ 
tion of roughly equal strength in two wavebands, in the 
optical/near-infrared and in the far-infrared at ~ 100 — 
200 //m. This equality implies that ~50% of the energy 
emitted by stars since the Big Bang has been absorbed by 
dust and is now in the far-infrared band. This is an av¬ 
erage value and studies of the relative cosmic evolution of 
the far-infrated and UV luminoisity density (Takeuchi, Buat 
and Burgarella 2005; Burgarella et al. 2013) imply that the 
fraction of the stellar light absorbed by dust is much higher 
at z > 1 during the epoch in which most of the stars in 
the Universe were formed. An important difference between 
the two peaks is that the sources that constitute the far- 
infrared background are rarer and consequently more lumi¬ 
nous. We show in this paper that in a small area of sky 
the optical/near-infrared background is composed of ~8000 
sources but the far-infrared background, which is of roughly 
equal strength, is composed of only ~1500 sources, imply¬ 
ing that the sources of the far-infrared background are a 
factor of ~5 more luminous than the more common sources 
composing the optical/near-infrared background. 

All of the HST studies above start from optical/near- 
infrared samples of galaxies. They all attempt to allow for 
the star-formation in each galaxy that is hidden by dust by 
using either the optical/near-infrared colours of the galaxy 


to estimate the extinction or, in two cases, the Spitzer 24- 
^m observations to estimate the far-infrared emission. The 
probem with the first method is that the luminous galaxies 
that make up the far-infrared background are so shrouded 
by dust that estimates of their star-formation rates from 
short-wavelength measurements are likely to be inaccurate. 
The problem with the second is the significant extrapolation 
necessary from the mid-infrared to the far-infrared. 

The second limitation of these studies is that all only 
reach conclusions about the morphologies of the galaxies in 
restricted samples, with the usual limits being redshift and 
stellar mass. As one example, Bruce et al. (2012) concluded 
that the star-forming galaxies in the redshift range 2 < 2 < 3 
and with stellar masses > 10 11 Mq have structures that are 
dominated by disks. Although conclusions like this are im¬ 
portant, they always leave open the possibility that most of 
the stars in the Universe were formed in a different range of 
redshift or in galaxies with stellar masses outside the lim¬ 
its of the sample. Our overall goal in this paper is to reach 
a more general conclusion about the fraction of the stellar 
mass-density that was formed in galaxies of different mor¬ 
phological types over the entire history of the Universe. 

To reach this conclusion, we have used a calorimetry 
technique, which is not new (Pagel et al. 1997; Lilly and 
Cowie 1987) but is particularly timely because of one of the 
successes of the Herschel Space Observatory. The method is 
based on the fact that the extragalactic background radia¬ 
tion is the repository of all the radiation emitted as the result 
of nuclear fusion in stars. For the method to be practical, it is 
necessary to resolve this background radiation into individ¬ 
ual sources. Deep surveys with the Hubble Space Telescope 
had already resolved the optical/near-infrared background 
into individual sources, and now the deepest Herschel sur¬ 
vey has also achieved this goal in the far infrared (Elbaz et 
al. 2011; Magnelli et al. 2013; Leiton et al. 2015). 

One of the advantages of the calorimetry technique is 
that it does not require us to estimate the star-formation 
rate in individual galaxies, which is difficult enough even 
for galaxies at zero redshift (Kennicutt 1998; Kennicutt 
and Evans 2012). For example, one method of estimating 
the star-formation rate in a galaxy is by assuming the far- 
infrared emission represents the bolometric luminosity of 
newly formed stars (Kennicutt 1998; Kennicutt and Evans 
2012). This is probably a reasonable assumption for the lu¬ 
minous galaxies that constitute the far-infrared background 
but it is clearly wrong in galaxies with low star-formation 
rates, in which much of the far-infrared emission is from dust 
heated by the old stellar population (Bendo et al. 2015). The 
calorimetry method avoids the need to make assumptions 
like this by linking the total energy emitted by a galaxy to 
the energy produced by nuclear fusion. 

Our secondary goal in this paper is to determine how 
the stellar mass-density in the Universe today depends on 
galaxy morphology. By comparing this distribution with the 
result from the caloriometry technique on the morphologies 
of the galaxies in which most stars were formed, it is possi¬ 
ble to quantify how important the transformation of galaxy 
morphology has been in the evolution of the galaxy popula¬ 
tion. 

We have achieved this secondary goal using two new 
wide-field surveys of the nearby Universe. The Herschel As- 
trophysical Terrahertz Large Area Survey (the Herschel- 
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ATLAS or H-ATLAS, Eales et al. 2010) is the Herschel 
extragalactic survey covering the largest area of sky, 550 
square degrees, roughly one eightieth of the entire sky. The 
survey consists of images at 100, 160, 250, 350 and 500 pm 
of five fields, two large fields at the North and South Celes¬ 
tial Poles (NGP and SGP) and three smaller fields on the 
celestial equator. The NGP and the equatorial fields were 
covered in the Sloan Digital Sky Survey, but more impor¬ 
tantly the equatorial fields have also been covered by a much 
deeper spectroscopic survey, the Galaxy and Mass Assembly 
(GAMA) redshift survey (Driver et al. 2011). The GAMA 
team has produced many different data products for the 
galaxies in their survey, and we have made liberal use of 
these data products in our investigation. 

One of the challenges in an investigation like this is to 
find a quantitative measure of galaxy morphology that can 
be used over a large range of redshift. The obvious way to 
measure the position of a galaxy on the morphological spec¬ 
trum from early-type to late-type is to measure the relative 
proportions of the stellar mass that is in the bulge/spheroid 
and the disk. However, there are two major practical prob¬ 
lems with doing this. The first is that although for many 
years it was believed that the spheroids of galaxies are de¬ 
scribed by the de Vauclouleur R 1,/4 model (de Vaucouleurs 
1948), it is now realised that galaxy bulges follow a range of 
intensity profiles (Graham 2013 and references therein). Sec¬ 
ond, estimates of the proportion of the galaxy’s stellar mass 
that is in the bulge depend strongly on the assumptions 
made about the dust extinction within the galaxy (Driver 
et al. 2008; Graham and Worley 2008), which requires great 
care and spatial resolution to measure accurately. 

We have instead used the simpler approach of using the 
Sersic index as a measure of a galaxy’s morphology. The 
Sersic function is 

I( r ) = I e exp(-b n [(-) 1/n -1}) (1) 

r e 

in which I e is the intensity at the effective radius, r e , that 
contains half of the total light from the model, n is the Sersic 
index and b„ is a function that depends on n (Sersic 1963; 
Graham and Driver 2005). This useful function can be used 
to fit the structures of all galaxy types. If the value of the 
Sersic index is one, the Sersic profile becomes an exponen¬ 
tial intensity profile, which in the low-redshift Universe is 
the intensity profile followed by a disk. If the value of the 
Sersic index is four, the Sersic profile becomes the de Vau¬ 
clouleur R 1 ' 4 profile, which is the intensity profile followed 
by ellipticals in the low-redshift Universe. The advantage of 
adopting this method of measuring the structure of a galaxy 
compared to measuring the bulge-to-disk ratio is that it is 
much simpler, not requiring any assumptions about the ef¬ 
fect of dust extinction or the bulge profile (Graham and 
Worley 2008). 

In Section 2, we describe our investigation of how the 
stellar mass-density and star-formation rate per unit comov¬ 
ing volume depend on galaxy morphology in the Universe 
today. Section 3 describes the caloriometry method that we 
have used to estimate how the mass of stars that has formed 
over the history of the Universe depends on the morphology 
of the galaxy in which those stars were formed. Section 4 
gives the results of the calorimetry method. The results and 


their implications are discussed in §5 and the conclusions are 
given in §6. We assume the cosmological parameters given 
from the Planck 2013 cosmological analysis (Planck Collabo¬ 
ration XVI 2014): a spatially-flat universe with Dm = 0.315 
and a Hubble constant of 67.3 km s _1 Mpc -1 . 


2 THE UNIVERSE TODAY 

2.1 The relationship between star-formation 
density and morphology 

In order to carry out a census of star formation, it seems sen¬ 
sible to start from a far-infrared survey, since newly formed 
stars are hidden by dust from optical telescopes. This proce¬ 
dure might lead to the star-formation rate being underesti¬ 
mated if there is a population of galaxies in which stars are 
forming but in which there is little dust. However, we show 
in an accompanying paper that at low redshifts we do not 
appear to be missing any star-forming but dust-free galaxies 
(Eales et al. 2015). 

We have therefore used as our starting point the far- 
infrared H-ATLAS survey rather than the optical GAMA 
survey. We have used the data for the H-ATLAS equato¬ 
rial fields (Eales et al. 2010), which are also covered by the 
GAMA survey and have a total area of ~ 160 deg 2 . The ver¬ 
sion of the H-ATLAS data used in this paper is the Phase 
1 version 3 internal data release, which was obtained us¬ 
ing very similar methods to those used for the H-ATLAS 
Science Demonstration Phase, rather than the public data 
release described by Valiante et al. (2015) and Bourne et 
al. (2015). The maps at 250, 350 and 500 pm were made 
using the methods described by Pascale et al. (2011) and 
the maps at 100 and 160 pm were made using the method 
described in Ibar et al. (2010). The detection of the sources 
was performed on the 250-pm images using the methods de¬ 
scribed in Rigby et al. (2011). Flux measurements at 350 and 
500 pm were generally made by measuring the fluxes at the 
250-pm positions on the 350-pm and 500-pm images after 
convolving the latter images with the point spread function 
of the telescope. We measured fluxes on the less sensitive 
100 and 160-pm images using aperture photometry, with an 
aperture diameter chosen to maximise the signal-to-noise for 
point sources. For the small percentage of H-ATLAS sources 
associated with nearby galaxies, we measured flux densities 
by carrying out aperture photometry at all wavelengths with 
an aperture diameter chosen to recover, as nearly as possi¬ 
ble, all the flux. 

The Phase 1 catalogue consists of the 109,231 sources 
detected at > 5<r at either 250, 350 or 500 pm (the maps at 
100 and 160 pm are much less sensitive than those at longer 
wavelengths). These signal-to-noise limits corresponds to 
flux densities at 250, 350 and 500 pm of 32, 36 and 45 mjy, 
respectively. 

As in the Science Demonstration Phase, we used a 
likelihood-ratio analysis (Smith et al. 2011) to match the 
sources to galaxies in the Sloan Digital Sky Survey (SDSS, 
DR7; Abazajian et al. 2009) with r < 22.4. With this 
method, we can estimate the probability (the reliability - 
R) of a galaxy being the true association of a submillimetre 
source. We treat galaxies with R > 0.8 as being likely coun¬ 
terparts to the submillimetre sources. There are counter¬ 
parts with R > 0.8 for 29,761 of the 104,657 sources detected 
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at 5cr at 250 /im. Of these counterparts 14,917 have spectro¬ 
scopic redshifts, mostly from the GAMA redshift survey but 
some from the SDSS. For the counterparts, there is matched- 
aperture photometry in u, g, r, i and z from the SDSS; in Y, 
J, H and K from the UKIRT Large Area Survey, which is 
part of the UKIRT Infrared Deep Sky Survey (Lawrence et 
al. 2007); and in near-ultraviolet and far-ultraviolet filters 
from the Galaxy Evolution Explorer (Martin et al. 2005). 
We have used the photometry to estimate photometric red- 
shifts using the artificial neural network code ANNZ (Smith 
et al. 2011). 

Our base sample for the analysis in this paper is the 
2250 Herschel sources detected at > 5 a at 250 /im in the 
fields common to H-ATLAS and GAMA that have reliable 
optical counterparts with spectroscopic redshifts < 0.1. The 
fraction of Herschel sources for which we can find counter¬ 
parts using the likelihood technique declines with the red- 
shift of the source, but at z < 0.1 the fraction of sources with 
identified counterparts should be close to 100% (Smith et al. 
2011). Virtually all the counterparts below this redshift have 
spectroscopic redshifts (Smith et al. 2011). 

We have estimated the key properties of each galaxy, 
in particular its star-formation rate, from the galaxy’s spec¬ 
tral energy distribution—the 11 photometric measurements 
listed above plus the five measurements in the Herschel 
bands—using the galaxy model MAGPHYS (Da Cunha, 
Chariot and Elbaz 2008). Some of the characteristics of 
MAGPHYS important for our analysis are that it takes ac¬ 
count of both the star formation that is seen directly in 
the optical waveband and the star formation that is hidden 
by dust, and that the model spectral energy distribution 
is fit to the observed spectral energy distribution using a 
probabilistic method, which explores a very large number 
of possible star-formation histories and leads to a probabil¬ 
ity distribution for each galaxy property. Further details of 
the application of MAGPHYS to the H-ATLAS galaxies are 
given in Smith et al. (2012b). We use the median estimate 
of the star-formation rate for each galaxy, SFRi, derived 
from the probability distribution of the star-formation rate 
produced by MAGPHYS. 

Our estimate of the star-formation rate per unit comov¬ 
ing volume in the Universe today is then 

q p jd 

SFR density = EE^ 1 ( 2 ) 

z 

in which Vi is the accessible volume in which the galaxy 
could have been detected in the HATLAS survey. In calcu¬ 
lating the accessible volume, we assume that the spectral 
energy distribution of the galaxy in the Herschel bands fol¬ 
lows a modified black body, with the dust temperature and 
emissivity index both a linear function of the logarithm of 
the 250-/mi monochromatic luminosity (Eales et al. in prepa¬ 
ration). However, at such low redshifts our estimate of the 
accessible volume of a galaxy depends only very weakly on 
the assumption made about its spectral energy distribution. 

The GAMA team has measured the value of the Sersic 
index for all these galaxies from images in nine different 
passbands (Kelvin et al. 2012). They find that there is a 
small increase in the measured value of the Sersic index when 
the wavelength of the image used to measure the Sersic index 



Sersic index 

Figure 1. Star-formation rate per unit comoving volume in the 
Universe today as a function of Sersic index. 

is increased. We have used the values they have measured 
from the SDSS r-band images, since this band corresponds 
well to the typical rest-frame wavelength of the images used 
to estimate the Sersic indices of the high-redshift galaxies 
(§3.5). We have measured how the star-formation rate per 
unit comoving volume in the Universe today depends on the 
Sersic index of the galaxy in which the stars are forming by 
restricting the sum in equation 2 to the galaxies that fall in 
each bin of Sersic index (Fig. 1). The figure shows that today 
the star-formation rate per unit comoving volume peaks at 
a Sersic index of ~1, implying that, today in the Universe, 
most stars form in galaxies that are dominated by a disk. 


2.2 The relation between stellar mass-density and 
morphology 

We can estimate the stellar mass-density as a function of 
Sersic index in a similar way. This time we have started 
from the optical GAMA survey, since there are many ETGs 
that contain large stellar masses but which are forming stars 
at a slow rate, have little dust, and are difficult to detect 
with Herschel. We have started from the GAMA I database, 
which includes galaxies with r-band magnitudes brighter 
than 19.4 in the GAMA9 and GAMA15 fields and 19.8 in the 
GAMA15 fields. We exclude galaxies with 2 < 0.002, since 
their velocities are likely to have a significant component 
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from peculiar motions. Our sample is the 15724 galaxies with 
spectroscopic redshifts in the redshift range 0.002 < z < 0.1. 
The stellar mass-density in the Universe today is given by 


Stellar mass density = 


Mti 


(3) 


in which M*i is the stellar mass of the i’th galaxy and Vi 
is again accessible volume. We have used the stellar masses 
for the galaxies derived by Taylor et al. (2011), since there 
are not yet MAGPHYS estimates of the stellar mass for all 
the GAMA galaxies. For the H-ATLAS galaxies, for which 
there are estimates from both methods, there is a system¬ 
atic offset of 0.17 in logio(stellar mass) between the two 
estimates, but once this offset is corrected for, the agree¬ 
ment is good, with a root-mean-squared discrepancy of 0.13 
in logio(stellar mass) between the two estimates (Eales et 
al. 2015). The offset is not important for this analysis in this 
section. We have calculated the accessible volume for each 
galaxy using an optical spectral energy distribution derived 
for each galaxy from its multi-band photometry (Taylor et 
al. 2011). We can estimate the dependence of stellar mass 
density on Sersic index by restricting the sum in equation 3 
to the galaxies that fall in each bin of Sersic index. Figure 
2 shows the stellar mass-density in the Universe today as a 
function of Sersic index. 

As we discussed in §1, there is no agreement about the 
best way to separate ETGs and LTGs. In this paper, we use 
the operational definition that ETGs and LTGs are galaxies 
with a Sersic index greater than and less than 2.5, respec¬ 
tively - a definition which has been used often in the liter¬ 
ature (e.g. Shen et al. 2003; Barden et al. 2005; Buitrago 
et al. 2013). One argument suggests that it is a definition 
which has some physical meaning. Figure 1 shows that the 
distribution of star-formation rate per unit comoving vol¬ 
ume is fairly symmetric, with a peak at n = 1. This value of 
n is the one expected if the structure of the galaxy is dom¬ 
inated by an exponential disk, suggesting that the physical 
meaning of the peak is that star formation mostly occurs 
in disk-dominated galaxies. Therefore it seems reasonable 
to devise an ETG/LTG classification that puts this peak in 
one class. With a dividing line of n = 2.5, the galaxies in 
this peak all fall in the class of LTGs, but with a lower value 
of n some of the galaxies in this peak would be classified 
as ETGs; and with a higher value of n some of the galaxies 
that are outside this peak would be classified as LTGs. Nev¬ 
ertheless, despite this argument, our definition is basically 
an operational one that allows us to compare the morpholo¬ 
gies of galaxies in the Universe today with those of galaxies 
in the distant past. Using this definition, we calculate from 
the distribution shown in Figure 1 that ~89% of the star 
formation in the Universe today is contained in LTGs. 

Using this definition, we calculate from Figure 2 that 
51% of the stellar mass-density today is in ETGs. Kelvin et 
al. (2014a) have used the optical images of galaxies in the 
GAMA survey to classify the galaxies into different Hubble 
types. Kelvin et al. (2014b) then used these classifications to 
show that 34% of the stellar mass-density today is in ellipti¬ 
cal galaxies and 37% is in SO-Sa galaxies, and so 71% of the 
stellar mass-density is in galaxies that would be tradition¬ 
ally assigned to the ETG class. The fact that our estimate 



Sersic index 

Figure 2. Stellar mass per unit comoving volume in the Universe 
today as a function of Sersic index. 


is lower is most likely explained by the fact that many of 
the Sa galaxies and and some SOs have structures that are 
dominated by disks (Laurikainen et al. 2010). 

Another way to classify galaxies is to separate them 
using their colours into star-forming galaxies and ‘passive’ 
or ‘quiescent’ galaxies in which star formation has largely 
stopped (e.g. Peng et al. 2010). Although we present argu¬ 
ments in our accompanying paper (Eales et al. 2015) for 
why this classification is misleading and hides some of the 
physics that is going on in the galaxy population, it is true 
that there is a large overlap between the red passive quies¬ 
cent galaxies and ETGs and between the blue star-forming 
galaxies and LTGs. Using the stellar mass functions for star¬ 
forming and passive galaxies given in Baldry et al. (2012), 
we calculate that 73% of the stellar mass-density today is in 
passive galaxies. The obvious explanation of our lower esti¬ 
mate of 51% for the fraction of the stellar mass-density that 
is contained in ETGs is if there is a population of passive 
galaxies which have structures that are dominated by disks 
(Van den Bergh 1976; Cortese 2012). 
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3 STAR FORMATION OVER THE HISTORY 
OF THE UNIVERSE - METHOD 

3.1 The calorimetry equation 

We make the assumption that all the extragalactic back¬ 
ground radiation is due to nuclear fusion in stars. A number 
of theoretical arguments and observational results show that 
~10% of the background radiation comes from the release 
of gravitational energy in active galactic nuclei (Dwek and 
Krennrich 2013), but this should have very little effect on 
our final results since the absolute value of the intensity of 
the background is not used in our method. 

The basic equation we use links the background radia¬ 
tion, B v , produced as the result of nuclear fusion in stars, 
with the elements that are produced by the nuclear fusion: 



0.007pc 3 u 
4tt(1 + z) 


(4) 


In this equation p = p(z ) is the average co-moving density of 
elements produced by nuclear fusion in stars per unit time 
at a redshift 2 , t is cosmic time, tH is the age of the Uni¬ 
verse and the numerical factor, 0.007, is the fraction of the 
mass of hydrogen that is turned into energy when hydrogen 
is fused into helium. Since the products of nuclear fusion 
accumulate over time, this equation allows us to relate the 
energy dumped into the background radiation by nuclear fu¬ 
sion at a redshift 2 to the products of this nuclear fusion in 
the Universe today. The calorimetry equation is exact and 
can be derived directly from the Robertson-Walker metric, 
and therefore does not require any assumptions about the 
values of the standard cosmological parameters. 

Although the equation is exact, for it to be a useful 
practical tool it is necessary to make some approximations. 
One that has been used in the past (Pagel 1997) is to as¬ 
sume that the density today of all the elements produced by 
nucleosynthesis, which includes helium, is proportional to 
the density today of the elements heavier than helium (the 
metals). In this paper, we have tried the different approach 
of connecting the stellar mass-density in the Universe to¬ 
day with the energy produced from nucleosynthesis in these 
stars that is now in the background radiation. 

With a few assumptions, which we will discuss below, 
we can write the following two equations: 


T h 

p* = (1 — R) / SFR{t)dt 
Jo 


(5) 


C t h 


roo r 00 

SFR(t)dt = 

Jo Jo 


4-7t(1 + z)J(z, u) 
0.007c 3 


dudz 


( 6 ) 


in which p* is the current average mass-density of stars in 

the Universe, SFR(t) is the mass of stars produced per unit 
volume per unit time at a redshift 2 , R is the fraction of 
that newly formed stellar mass that is subsequently ejected 
into the interstellar medium as the result of stellar winds 
(e.g. Leitner and Kravtsov 2011), J[z , u) is the background 
radiation produced by objects at redshift 2 in a unit interval 
of redshift, and th is again the current age of the Universe. 


A basic assumption that we have made to produce the 
first equation is that the stars formed in the Universe at 
a time t are still in the Universe today. This assumption 
seems a reasonable one since most of the mass in a stellar 
population is in long-lived low-mass stars. 

An assumption made to produce the second equality is 
that all the available nuclear energy in stars has been lib¬ 
erated and dumped into the background radiation by the 
current epoch, which is obviously not true because the Sun 
is still shining. However, it is not a crazy assumption be¬ 
cause most of the energy produced by a population of stars 
is produced by the high-mass stars, which do burn through 
their available nuclear energy very quickly. This is confirmed 
by Figure 3, where we have used models of stellar popula¬ 
tions (Maraston 2006) to predict how the total energy that 
has been emitted by a single stellar population depends on 
age for various assumptions about the initial mass func¬ 
tion (IMF). The figure shows that for a Salpeter or Kroupa 
IMF most of the energy has been emitted by 1 Gyr and the 
amount of energy that has been radiated (and dumped into 
the background radiation) increases by only ~50% as the 
population ages from 1 Gyr to 10 Gyr. 

Our main results in this paper, however, were not ob¬ 
tained by assuming either equality but by combining the two 
equations in the form of a proportionality in which we as¬ 
sume that the current stellar mass-density in a class of galax¬ 
ies, p* c , is proportional to the contribution of the galaxies 
in that class to the background radiation: 


p*c oc 



4-7t(1 + z)J(z , u) 
0.007c 3 


dudz 


(7) 


The assumptions behind this proportionality are slightly dif¬ 
ferent from the assumptions behind the equalities. We are 
now assuming that the fraction of newly formed stellar mass 
lost by stellar winds and that the fraction of the available 
nuclear energy liberated by the current epoch (see above) 
do not vary between galaxy classes. The second quantity is 
probably different between ETGs and LTGs, but we argue 
in §5 that the difference actually strengthens the main result 
in this paper. 

One further assumption that is built-in to the calori¬ 
metric method is that the IMF does not vary significantly 
between galaxy classes. In Figure 3, we have explored the 
importance of this assumption by predicting the relation¬ 
ship between bolometric luminosity and time for different 
assumptions about the IMF. The figure shows that the dif¬ 
ference between a Kroupa and a Salpeter IMF is quite small, 
but our results would be changed significantly if some galaxy 
class had a bottom-heavy IMF, because these galaxies might 
then contain a large mass of stars today without having 
contributed much radiation to the extragalactic background 
radiation. 


3.2 The extragalactic background radiation 

The extragalactic background radiation (Dole et al. 2006; 
Dwek and Krennrich 2013) is dominated by the radia¬ 
tion in two wavebands of roughly equal strength: in the 
optical/near-infrared, where the background peaks at ~1- 
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log age/yr 


Figure 3. The total energy that has been emitted by a given 
age for a single stellar population with a mass of 10 11 Mq. We 
have used the models of stellar populations of Maraston (2006). 
The different lines are for different assumptions about the stellar 
initial mass function. ‘Kroupa’ and ‘Salpeter’ refer to initial mass 
functions measured in the solar neighborhood, whereas ‘top-heavy 
and ‘bottom-heavy refer to initial mass functions with a larger 
proportion of high-mass and low-mass stars, respectively. 


2 pm , and in the far infrared, where the background peaks 
at ~ 100-200 /im. 

The far-infrared background radiation has been mea¬ 
sured by the FIRAS and DIRBE experiments on the Cosmic 
Background Explorer satellite (COBE) (Fixsen et al. 1998; 
Dole et al. 2006; Fig. 4). In calculating the contribution of 
individual sources to the background (§3.3), at 160, 250, 350 
and 500 pm we used the FIRAS observations of the back¬ 
ground. The l-<7 uncertainty on the FIRAS estimate of the 
background is ~30%. At 100 pm we used the DIRBE mea¬ 
surement placed on the FIR AS flux scale (Dole et al. 2006), 
since the FIRAS measurements did not extend down to this 
wavelength. 

In the optical/near-infrared, there have been three dif¬ 
ferent techniques used to estimate the strength of the extra- 
galactic background radiation: 1) estimates from integrat¬ 
ing the galaxy counts; (2) absolute measurements of the 
background using the DIRBE instrument on COBE (Cam- 
bresy et al. 2001; Levenson et al. 2007); (3) searches for 
absorption effects in the gamma-ray spectra of distant ac¬ 
tive galactic nuclei (Dwek and Krennrich 2013). The second 
method would be ideal except for the systematic effects on 
the DIRBE near-infrared measurements from sunlight scat- 
ted by interplanetary dust. Even after a correction has been 
made for this scattered light, there is generally a clear rela¬ 
tionship between the measurements of the background radi¬ 
ation and ecliptic latitude (Cambresy et al. 2001; Levenson 
et al. 2007), implying a residual contamination from inter¬ 
planetary dust (reliable determinations will probably require 
measurements from beyond 3 A.U. from the Sun, where the 


contamination from scatted light is small (Levenson et al. 
2007)). These absolute measurements are 2-4 times greater 
than the estimate from integrating the galaxy counts, but 
upper limits from the gamma-ray spectra of active nuclei are 
generally lower than these absolute measurements, also sug¬ 
gesting that the DIRBE measurements are too high (Dwek 
and Krennrich 2013). For these reasons we have taken as 
our estimate of the background at 1.6 pm the integral of the 
deep 1.6-pm ( Hab ) galaxy counts made with Wide Field 
Camera 3 on the Hubble Space Telescope (Windhorst et al. 
2011 ). 

3.3 Samples of the background radiation 

To apply the calorimetry technique (§3.1), we ideally need 
to be able to resolve the extragalactic background radiation 
into individual sources. By measuring the structures of the 
galaxies producing the emission, we can then use equation 
7 to determine the proportion of stars that have formed in 
galaxies with different structures. 

In the southern field of the Great Observatories Ori¬ 
gins Deep Survey (GOODS-South), we can get surprisingly 
close to this ideal. This field has been surveyed by the deep¬ 
est survey with the Herschel Space Observatory, Herschel- 
GOODS (Elbaz et al. 2011; Magnelli et al. 2013) at 100 and 
160 pm, close to one of the two peaks in the extragalac¬ 
tic background radiation. It has also been observed as part 
of the Cosmic Assembly Near-Infrared Deep Extragalactic 
Legacy Survey (CANDELS), an optical/near-infrared sur¬ 
vey carried out with the Hubble Space Telescope with a 
main wavelength of 1.6 pm, close to the other peak of the 
extragalactic background radiation (Koekemoer et al. 2011; 
Guo et al. 2013). The CANDELS team has measured Sersic 
indices for the galaxies detected in the survey (Van der Wei 
et al. 2012) from data that are well matched to the ground- 
based images used to estimate the Sersic indices in our inves¬ 
tigation of the Universe today (§2, §3.5). We now show how 
these surveys can be used to resolve the two peaks in the 
extragalactic background radiation into individual sources. 
Although the total area covered by the surveys of GOODS- 
South is only ~100 arcmin 2 , we show in §4 that cosmic vari¬ 
ance is not a significant problem, and thus that the samples 
of sources provided by these surveys constitute fair samples 
of the extragalactic background radiation. 

3.3.1 Sampling the far-infrared background with the 
sample of Herschel sources 

The public catalogue for Hcrschel-GOODS (H-GOODS) is 
available at http://hedam.lam.fr/GOODS-Herschel/ This 
catalogue was obtained by starting with the assumption that 
the sources in the 100 and 160-pm images are at a set of 
‘prior positions’ drawn from a catalogue of sources detected 
at 24 pm with the Spitzer telescope, and then varying the 
100-pm and 160-pm flux densities of the sources to match 
the observed structure in the images (Elbaz et al. 2011; 
Magnelli et al. 2013). The Spitzer catalogue itself was ob¬ 
tained in a similar way from a set of prior positions drawn 
from a catalogue of sources detected at 3.5 pm with the 
Spitzer telescope (Elbaz et al. 2011; Magnelli et al. 2013). 
This technique has the advantage that it produces very accu¬ 
rate positions for the H-GOODS sources because the 3.5-pm 
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sources have much more accurate positions than is possible 
to obtain with Herschel observations by themselves. How¬ 
ever, it has the potential disadvantage that it could miss 
Herschel sources if they are not present in the prior cat¬ 
alogue. By comparing the H-GOODS catalogue with one 
obtained without using the prior catalogue, the H-GOODS 
team estimates that the catalogue made using the prior po¬ 
sitions misses <4% of the sources at 160 pm (Magnelli et al. 
2013). Furthermore, any missing sources will be among the 
fainter ones, and thus should have a negligible effect on the 
calorimetric technique. 

The H-GOODS catalogue contains sources that were 
detected at 100 and 160 /im at > 3cr. Sample 1 in our anal¬ 
ysis contains the 527 sources detected at > 3a at 160 pm 
in the region of GOODS-South also covered by the CAN- 
DELS survey. To estimate the contribution of this sample 
to the background radiation, we need to make measurements 
of the flux density at 100, 160, 250, 350 and 500 /im for ev¬ 
ery source in the sample; even if an individual detection is 
below the 3 a level, the result of summing the results for all 
the sources is often highly significant. We essentially carry 
out a statistical ‘stacking analysis’ to estimate the contribu¬ 
tion from all the sources, even those that are not detected 
individually at a significant level (Marsden et al. 2009). 

There are no H-GOODS images at the three longer 
wavelengths, so at these wavelengths we measured the flux 
densities of the sources from images obtained as part of the 
HERMES survey (Oliver et al. 2012). To measure the flux 
density of a source, the standard technique is to first con¬ 
volve the image with the point spread function of the tele¬ 
scope, since this maximizes the signal-to-noise of the source. 
However, this is only correct if the noise on the image is un¬ 
correlated between pixels, which is true of the instrumental 
noise on the HERMES images. Most of the noise on the 
HERMES images, however, is confusion noise, the noise due 
to faint submillimetre sources that are too faint to detect 
individually, and this noise is correlated between pixels. We 
did not therefore convolve the HERMES images but simply 
measured the flux densities from the pixel in each HERMES 
image that contained the position of the H-GOODS source. 

The background level on a Herschel image is usually 
close to but not exactly zero. To estimate the residual back¬ 
ground emission and the errors on the flux densities, we used 
the Monte-Carlo approach of measuring flux densities at a 
large number of random positions on the HERMES images 
without making a correction for any residual background 
emission. We found background levels of -0.48±0.08 mjy, - 
0.74±0.08 mjy and -1.1±0.08 mjy at 250, 350 and 500 /im, 
respectively, which we used to correct our measurements of 
the flux densities of the H-GOODS sources. From the Monte- 
Carlo simulation, the errors on these flux densities at 250, 
350 and 500 pm are 6.2 mjy, 6.5 mjy and 6.2 mjy, respec¬ 
tively, very similar to the Herschel confusion noise at these 
wavelengths (Nguyen et al. 2010). 

We used the 160-pm flux densities from the catalogue, 
but there was often not a 100-pm measurement in the cat¬ 
alogue. If a measurement did not exist in the catalogue, 
we used aperture photometry to measure the flux den¬ 
sity from the H-GOODS 100-pm image, using an aperture 
with a radius equal to the full-width-half-maximum of the 
point spread function (6.7 arcsec). We again used a Monte- 
Carlo approach to estimate the residual background level 


on the image and the errors on the flux densities, perform¬ 
ing aperture photometry at a large number of random po¬ 
sitions on the image without correcting for a residual back¬ 
ground. We found a residual background level in our aper¬ 
ture of 0.3±0.03 mjy, which we used to correct the aperture 
photometry of the H-GOODS sources. For the H-GOODS 
sources for which 100-pm photometry exists in the public 
catalogue, we compared our measured flux densities with 
the catalogued flux densities to derive an aperture correc¬ 
tion to bring our measurements on to the same flux scale as 
the catalogue. As a result of this comparison, we have mul¬ 
tiplied our aperture flux densities by a factor of 1.82. From 
the Monte-Carlo simulation, the error on our measured flux 
densities is 2.4 mjy. 

We estimated the contribution of the sample to the 
background radiation at each wavelength by adding the flux 
densities of the individual sources. We estimated the error 
on the sum by adding in quadrature the errors on the in¬ 
dividual flux densities. The errors on the total flux density 
of the sources are 2%, 0.6%, 4%, 5% and 10% at 100, 160, 
250, 350 and 500 pm, respectively, much less than the er¬ 
rors on the measurements of the background radiation at 
these wavelengths. Thus the errors on the percentage of the 
background radiation resolved by the sources are dominated 
by the errors on the background measurements themselves, 
which we assume are 30% for the measurements with FIRAS 
at the four longest wavelengths and 44% for the DIRBE mea¬ 
surement at 100 pm (Dole et al. 2006). Using these errors, 
our estimate of the fraction of the background contributed 
by the H-GOODS sample is 0.551 q ^ at 100 pm, 0.67+Qig 
at 160 pm, 0.481 q'ij at 250 pm, OAOIq'^ 350 pm, and 
0.451 q'iq at 500 pm (Fig. 4). The H-GOODS team has es¬ 
timated the fraction of the background resolved by the H- 
GOODS sources as 75% at both 100 and 160 pm (Magnelli 
et al. 2013). There is therefore a sizeable difference at 100 
pm for which we do not know the reason. 

The sample of sources detected with Herschel in the H- 
GOODS survey therefore constitutes a significant fraction of 
the background radiation but still falls significantly short, 
whether we use our estimates or the estimates of the H- 
GOODS team, from resolving all the background radiation. 

3.3.2 Sampling the far-infrared background with a sample 
of Spitzer sources 

We show now that we can actually sample a larger fraction 
of the far-infrared background using a sample of sources de¬ 
tected at 24 pm in the GOODS-South region by the Spitzer 
Space Telescope. This sample (sample 2) consists of 1557 
sources detected at 24 pm in the same area of GOODS- 
South that is covered by CANDELS. We measured the flux 
densities of each of these sources at 250, 350 and 500 pm 
from the HERMES images in the same way as for sample 1. 
At 100 and 160 pm, we used for preference the flux densities 
given in the H-GOODS public catalogue. If there was not a 
100-pm measurement in the public catalogue, we carried out 
aperture photometry on the H-GOODS 100-pm image in ex¬ 
actly the same way as for sample 1. For sample 2 this was 
often also necessary at 160 pm. In this case we measured 
a 160-pm flux density from the H-GOODS 160-pm image 
through an aperture with a radius equal to the full-width- 
half-maximum of the point spread function (11.0 arcsec). We 
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used the Monte-Carlo method of performing aperture pho¬ 
tometry at a large number of random positions on the H- 
GOODS 160-^im image to estimate the residual background 
level in this aperture of 2.0±0.07 mjy, which we allowed for 
in the measurements of the flux densities of the H-GOODS 
sources. By comparing our measured flux densities with the 
ones in the catalogue, where they exist, we derived a cor¬ 
rection factor of 2.0, by which we multiplied our aperture 
flux densities to bring them on to the same flux scale as the 
catalogue. From the Monte-Carlo simulation, our estimate 
of the error in our 160-pm flux densities is 5.6 mjy. 

For sample 2, the error in the total flux density of the 
sources is 5%, 7%, 4%, 5% and 8% at 100, 160, 250, 350 
and 500 pm, respectively, which are again much less than 
the errors in the level of the backround. Using 30% as the 
error in the FIRAS measurements of the background at the 
four longest wavelengths and 44% as the error in the DIRBE 
background measurement at 100 pm, we calculate that the 
sources in sample 2 compose 0.671°' jq of the background 
radiation at 100 pm, 0.901 q' 21 at 160 pm, 0.801 qi 9 at 250 
pm, 0.891 q'2i at 350 pm, and 0.901 q'2i at 500 pm. Given the 
uncertainties, the combined far-infrared flux density from 
sample 2 is fully consistent with the far-infrared background 
radiation over this range of wavelengths. The sources in this 
sample therefore represent a fair sample of the extragalactic 
background radiation in one of its two main peaks. 


3.3.3 Sampling the optical/near-infrared background with 
the HST 

The other peak in the extragalactic background radiation 
is in the optical/near-infrared waveband. As our sample of 
this peak (sample 3) we use the 8488 galaxies detected by 
CANDELS with Hab < 24.5 in the GOODS-South field. We 
have used this photometric band because its central wave¬ 
length (~1.6 pm) is close to the peak of the optical/near- 
infrared background. We have used the magnitude limit be¬ 
cause galaxies fainter than this do not generally have reli¬ 
able measurements of the Sersic index (Van der Wei et al. 
2012). This sample constitutes ~94% of the extragalactic 
background radiation at 1.6 pm. The high value is the re¬ 
sult of our assumption that the extragalactic background 
radiation at this wavelength is given by the integral of the 
galaxy counts (§3.2); since the slope of the galaxy counts 
flattens at Hab — 21 (Windhorst et al. 2011), the frac¬ 
tion of the background radiation produced by galaxies with 
Hab < 24.5 is very high. As long as there is not a substan¬ 
tial contribution to the extragalactic background radiation 
from starlight that is not detected in even the deepest expo¬ 
sures with the HST (§5), sample 3 constitutes a fair sample 
of the extragalactic background radiation in its second peak. 

We have included all the data for samples 1 and 2 that 
was used in this paper in on-line tables. The tables includes 
the positions of the Herschel and Spitzer sources, the flux 
densities at the five Herschel wavelengths, the spectroscopic 
and/or photometric redshifts of the sources, the positions of 
the counterparts on the HST images (§3.4) and the Sersic 
indices of the counterparts (§3.5). 



Wavelength//^, m 

Figure 4. The extragalactic far-infrared background radiation 
plotted against wavelength, and our estimates of the contribu¬ 
tion to this background from the samples described in this paper. 
The heavy solid line shows the measurement of the background 
with the FIRAS experiment on COBE (Fixsen et al. 1998), with 
the light dashed lines showing the ±lcr errors on the background. 
The black points show the measurements made by the DIRBE 
experiment on COBE (Dole et al. 2006). The red points show 
our measurements of the contribution to the background radi¬ 
ation from the sources individually detected at 160 gm by the 
H-GOODS survey (sample 1). The blue points show our mea¬ 
surements, using the Herschel images of GOODS South, of the 
contribution to the background radiation from the sources that 
were originally detected at 24 gm with Spitzer (sample 2). 


3.4 Identifying the counterparts to the Herschel 
and Spitzer sources 

There are measurements of the Sersic index for all the galax¬ 
ies in sample 3 by the CANDELS team (Van der Wei et al. 
2012). However, for the other two samples, before we can 
consider the structures of the galaxies producing the far- 
infrared emission, we need to identify these galaxies. 

The positions of the sources in samples 1 and 2 are ul¬ 
timately based on observations with the Spitzer Space Tele¬ 
scope at 3.5 /im, which have an angular resolution (FWHM) 
of 1.7 arcsec. The accuracy of the positions should there¬ 
fore be <1 arcsec, and we have made the assumption that 
any genuine counterpart to a Herschel/Spitzer source on the 
CANDELS H-band (1.6 gm) image must lie within 1.5 arc¬ 
sec of the position of the source. 

To assess the probability of a possible counterpart on 
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the CANDELS image within 1.5 arcsec of the position of 
the Herschel/Spitzer source being there by chance, we fol¬ 
low the frequentist technique commonly used in far-infrared 
astronomy of calculating a statistic, S, for the counterpart, 
defined as the number of CANDELS galaxies brighter than 
the 1.6-^tm flux density of the potential counterpart that 
are expected by chance to lie closer to the Herschel/Spitzer 
position than the potential counterpart (Dye et al. 2009). If 
there was more than one possible counterpart within 1.5 arc- 
sec, we chose the one with the lowest value of S as the most 
likely counterpart. We estimated the probability, P, of ob¬ 
taining such a low value of S, given the null hypothesis that 
there is no relation between the CANDELS galaxies and the 
Herschel/Spitzer sources, by the Monte-Carlo approach of 
measuring S for a large number of random positions on the 
image. Sample 1 contains 527 sources and we found possi¬ 
ble counterparts within 1.5 arcsec for 526. We estimated the 
number of counterparts that are there by chance by sum¬ 
ming P for the 526 counterparts, which is 3.0. Therefore 
<1% of the counterparts are likely to be spurious. Sample 
2 contains 1557 sources and we found possible counterparts 
within 1.5 arcsec for 1537. In this case, the sum of P is 9.9, 
again showing that <1% of the counterparts are likely to be 
spurious. 

From a formal statistical point of view, there is no rea¬ 
son why this frequentist technique should produce any bi¬ 
ases. However, one might worry that, because the flux at 
1.6 pm is used in the calculation, the technique is biased 
against highly obscured dusty galaxies, the kind that one 
might expect to be far-infrared sources. However, an inves¬ 
tigation shows that this is not the case. In sample 1, of the 
527 160-pm sources, there were 92 cases where there was 
more than one potential HST counterpart within 1.5 arc¬ 
sec. In only 5 of the 92 cases was there an object fainter 
than the selected counterpart but closer to the Herschel po¬ 
sition. A final argument that the method is reliable is that 
77% of the selected counterparts lie within 0.2 arcsec of the 
position of the Herschel source, whereas our Monte-Carlo 
analysis shows that the distribution of offsets for random 
positions reaches a maximum at >1 arcsec. 

3.5 The Sersic indices of the CANDELS galaxies 

In this section we consider whether any systematic errors 
might arise from comparing the Sersic indices of the H- 
ATLAS galaxies (§2) with those measured for the galaxies 
detected in CANDELS. 

The Sersic indices for the H-ATLAS sources were mea¬ 
sured from the r-band (0.62 pm) images from the Sloan 
Digital Sky Survey (SDSS) (Kelvin et al. 2012). The Sersic 
indices for the H-GOODS galaxies were measured from H- 
band (1.6 pm) images from the CANDELS HST programme 
(Van der Wei et al. 2012). The median redshift of the H- 
ATLAS galaxies is 0.07, whereas most of the H-GOODS 
galaxies lie at a > 1 (§4). For these redshifts, the HST im¬ 
ages of the H-GOODS galaxies (angular resolution of 0.17 
arcsec) and the SDSS images of the H-ATLAS galaxies (an¬ 
gular resolution ~1 arcsec) give almost exactly the same 
physical resolution, ~1.3 kpc. The r-band image of a galaxy 
at z = 0 and the H-band image of an H-GOODS galaxy 
at 2 = 1.6, which is fairly representative of the H-GOODS 
redshift distribution (§4), correspond to the same rest-frame 


wavelength. The similarity of the rest-frame wavelengths en¬ 
sures that there should be no biases created by the small 
changes that are evident when Sersic indices are measured 
at different wavelengths (Kelvin et al. 2012). The one slight 
mismatch is in the sensitivity to surface brightness for the 
two datasets. The HST image of an H-GOODS galaxy at 
2 = 1.6 has ~1 mag less surface-brightness sensitivity than 
the SDSS image at the same rest-frame wavelength of a 
galaxy at 2 = 0, and the difference increases at z > 1.6. 
Van der Wei et al. (2012) have investigated the effect of the 
signal-to-noise on the CANDELS images on the measured 
value of the Sersic index, concluding that there is little bias 
down to Hab = 23.0 but that by Hab = 24.0 galaxies with 
high values of the Sersic index (n > 3) will have measured 
values that are lower than the true values by 25% (see their 
Table 3). Therefore, some of the fainter objects in our high- 
redshift sample may have measured Sersic indices that are 
too low. 

3.6 Practical application of the calorimetry 
technique 

Subject to the assumptions discussed in §3.1, equation 7 
shows that the stellar mass-density in the Universe today 
associated with the star-formation in the i’th source in a 
fair sample of the background radiation is proportional to 
Ji(l + Zi), in which Zi is the redshift of the galaxy and Ji 
is the contribution of this galaxy to the extragalactic back¬ 
ground radiation Q. Therefore, given estimates of Ji and 2 ; 
for all the galaxies in a fair sample of the background ra¬ 
diation, it is straightforward to estimate the stellar mass- 
density in the Universe today produced as the result of the 
formation of stars over the history of the Universe. Given 
measurements of the Sersic index for each galaxy in a sam¬ 
ple, it is also straightforward to estimate how the stellar 
mass-density today depends on the Sersic index of the galaxy 
in which those stars were formed. It is important to note that 
these stars may well now be in a galaxy very unlike that in 
which the stars were formed, either as the result of the merg¬ 
ing of galaxies or some other process that has transformed 
the structure of the galaxy. 

Since the flux densities of individual galaxies are often 
detected with very low signal-to-noise or are even negative, 
we calculated Ji in both the optical/near-infrared and the 
far infrared using the following formula, which provides a 
robust estimate in this situation: 

Ji = f ( FijAi/j (8) 

j 

In this formula Fij is the flux density of the i’th galaxy 
in the j’th band and A Vj is the frequency range associated 
with that band. We adopted the very simple approach of as¬ 
suming that the boundaries between two bands, for example 

1 There is no causal association, of course. We are using the sam¬ 
ples to represent the energy generated by nucleosynthesis over the 
history of the Universe, but the starlight we are detecting in the 
early Universe has no causal connection with the stars we see to¬ 
day. We are implicitly assuming the cosmological principle that 
any region of the Universe looks like any other. 
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the Herschel bands centred at 100 and 160 fim, lie equidis¬ 
tant in log(frequency) between the two central wavelengths. 
For samples 1 and 2, we calculated Ji from the measured 
flux densities of each galaxy at 100, 160, 250, 350 and 500 
/im. For sample 3, we calculated J) from measurements of 
the flux density in the following photometric bands (Guo et 
al. 2013): U-band, ACS F606W, ACS F814W, WFC125W, 
WFC160W, IRACBl and IRACB2. 

For redshifts we used, in order of preference, spectro¬ 
scopic redshifts, the photometric redshifts that form part of 
the H-GOODS public data release and the catalogue of pho¬ 
tometric redshifts for the CANDELS galaxies (Li-Ting et al. 
2014). 

After losing ~1% of the sources for which we could not 
find a counterpart on the CANDELS image, samples 1 and 
2, contain 526 and 1537 galaxies, respectively. We removed 
~10% of the galaxies in each sample for which the flag for 
the Sersic measurement (Van der Wei et al. 2012) implies the 
measurement is not reliable, leaving 467 and 1358 galaxies 
in samples 1 and 2, respectively. Of the galaxies that are 
left, 352 (75%) in sample 1 and 903 (66%) in sample 2 have 
spectroscopic redshifts; for the remaining galaxies we used 
the photometric redshifts. 

In the optical/near-infrared, sample 3 contains 8488 ob¬ 
jects with Hab < 24.5. The small size of the field means 
that there is the potential for effects due to cosmic variance 
(§4). These effects will be largest at low redshifts, and in¬ 
deed there are a large number of galaxies at 2 ~ 0.05. To 
minimize these effects we excluded 123 galaxies with red¬ 
shifts < 0.06. We also excluded 15 very bright objects with 
Sersic indices of < 0.6, which are probably stars or quasars. 
Of the remaining objects, we eliminated ~20% for which the 
flag implies there is not a reliable measurement of the Ser¬ 
sic index (Van der Wei et al. 2012). Of the remaining 6753 
galaxies, 2482 (37%) have spectroscopic redshifts, and for 
the others we used the photometric redshifts. 


4 RESULTS 

Figure 5 shows the redshift distributions of the three sam¬ 
ples. Since GOODS-South is such a small field, an obvious 
concern is whether it is a fair sample of the Universe and, in 
particular, of the extragalactic background radiation. We es¬ 
timated the cosmic variance for the three samples using the 
on-line calculator described by Treventi and Stiavelli (2008). 
Using a mean redshift of 1.0 and a redshift interval of 2.0, 
we estimate that cosmic variance leads to errors in the total 
numbers of sources in samples 1, 2 and 3 of 7, 6 and 5%, re¬ 
spectively. Therefore, the samples as a whole are likely to be 
fair samples of the overall background radiation. When we 
apply the on-line calculator to the individual redshift bins, 
which have a width of 0.2 in redshift, we find that the uncer¬ 
tainty in the number of sources in each redshift bin is ~30% 
for samples 1 and 2 (left-hand panel of Figure 5). Hence the 
lumpy nature of the redshift distributions in Fig. 5 is prob¬ 
ably the result of large-scale structure. Note that there are 
well-known over-densities in GOODS-South at z = 0.66 and 
a = 0.735 (Adami et al. 2005). 

Using the values of Ji( 1 + Zi) and equation 6, we have 
estimated the star-formation history of the Universe: the 
star-formation rate per unit comoving volume as a func¬ 


tion of redshift. Figure 6 shows the star-formation history 
derived from the data for the three samples. Not surpris¬ 
ingly, at a given redshift the star-formation rate for sample 
2 is higher than for sample 1 because the sources in sam¬ 
ple 2 comprise a larger fraction of the background radiation. 
The star-formation histories derived from all three samples 
are rather noisy, the consequence of large-scale structure, 
but all three star-formation histories peak somewhere in the 
redshift range 1 < 2 < 3, consistent with the results of previ¬ 
ous attempts to derive the Universe’s star-formation history 
(e.g. Hopkins and Beacom 2006). 

Figure 7 shows the estimates from the three samples 
for the relationship between the stellar mass-density in the 
Universe today and the Sersic index of the galaxy in which 
those stars contributing to that mass-density were formed. 
We showed in §2 that in the Universe today 51% of the stellar 
mass density is in ETGs, which we empirically defined as 
galaxies with a Sersic index > 2.5. If the structures of these 
galaxies have not changed between the time the stars now in 
them were formed and the current epoch, we would expect 
Figure 7 to look very like Figure 2. On the contrary, the 
distributions for all three samples are strikingly similar to 
Figure 1, which showed how the star-formation rate in the 
Universe today depends on the Sersic index of the galaxy in 
which the those stars are being formed. Figure 1 showed, not 
unexpectedly, that, today, most stars form in LTGs. Figure 
7 shows that this has always been true, that even though 
most stars today are in ETGs, most stars formed in LTGs. 

We can quantify this very easily by measuring the pro¬ 
portions of the distributions in Figure 7 that lie on either 
side of the n = 2.5 dividing line. The result from both the 
samples of the far-infrared background, samples 1 and 2, is 
that ~87% of the stellar mass that is in the Universe to¬ 
day must have formed in LTGs. The result from the sample 
of the optical/near-infrared background, sample 3, is that 
~79% of the stellar mass today was formed in LTGs. To as¬ 
sess the importance of the measurement errors for the Ser¬ 
sic indices on these estimates, we carried out a bootstrap 
analysis for the smallest sample, sample 1. We generated 
1000 artificial versions of sample 1, using a random-number 
generator to produce from the measured Sersic index and 
measurement error for each galaxy a new ‘measurement’ in 
each sample for that galaxy. From these 1000 samples, we 
estimate the statistical uncertainty in our estimate of 0.87 
for the fraction of the stellar mass-density formed in LTGs 
as only 1.9 x 1CU 6 . We discuss other more important sources 
of error in the following section. 

Taking an average of the results for the two peaks in the 
background radiation, which is reasonable since they have 
roughly equal strength, our results imply that ~83% of the 
stellar mass today was formed in LTGs. This is much higher 
than our estimate of 49% for the percentage of the stellar 
mass-density that is today in LTGs. The big difference in 
the values implies that the morphological transformation of 
galaxies must be an important part of galaxy evolution. 


5 DISCUSSION 

The importance of the calorimetric approach is that it pro¬ 
duces a quantitative measurement of the stellar mass-density 
that has formed in galaxies with different morphologies 
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Figure 5. Redshift distributions of the samples of the far-infrared background radiation (left) and the optical/near-infrared background 
radiation (right). The red line in the left-hand panel shows the redshift distribution for the sample of sources detected by Herschel 
(sample 1) and the blue dashed line shows the redshift distribution for the sample of sources detected with Spitzer (sample 2). 


which covers all cosmic epochs and galaxy masses, since the 
extragalactic background radiation contains all the energy 
ever emitted as the result of nuclear fusion in stars. The 
calorimetric technique, if there are no practical problems— 
we will discuss some of the possibilities next—makes it im¬ 
possible for star formation to hide anywhere. 

One fundamental problem we can do nothing about is 
that the method is based on the assumption of a universal 
initial mass function. As we showed in §3.1, a population 
of stars that are formed with a bottom-heavy initial mass 
function may lead to a large stellar mass-density in the Uni¬ 
verse today but not have dumped much radiation into the 
extragalactic background radiation. 

A less fundamental problem, but one we can still do 
nothing about, is that the uncertainties on the absolute val¬ 
ues of the extragalactic background radiation still leave some 
places for star formation to hide; if the far-infrared back¬ 
ground is at the upper end of the range set by the errors 
(Fig. 4), or if the optical/near-infrared background has a 
component that is not detected in the deep HST surveys, 
our samples may comprise substantially smaller percentages 
of the background radiation than we have estimated. 

Third, we have missed some of the energy produced by 
nucleosynthesis in stars because we have only considered the 


background radiation over two limited ranges of wavelength 
around each of the two peaks in the background radiation 
(§4). This will only be a problem if the galaxies comprising 
the missing part of the background radiation have system¬ 
atically different morphologies from the galaxies in our sam¬ 
ples of the two peaks. This is possible, but seems unlikely 
because the two peaks are separated by a factor of ~100 in 
wavelength, yet the morphological makeup of the galaxies 
in the samples of the two peaks are quite similar (compare 
the right and left panels of Fig. 7). 

Fourth, there is the issue that for any population of 
stars, not all the available nuclear energy will have been 
dumped into the extragalactic background radiation by the 
current epoch (§3.1, §4). ETGs generally have older stellar 
populations than LTGs (Kennicutt 1998) and thus should 
have dumped a larger fraction of the total available nuclear 
energy into the background radiation by the current epoch, 
and so if there is a systematic bias it will be in the sense 
that samples of the extragalactic background radiation pref¬ 
erentially contain ETGs. Therefore, this may have led us to 
underestimate the true difference beween the fraction of the 
stellar mass-density that is in LTGs today and the fraction 
of the total stellar mass-density that was formed in LTGs. 

Fifth, there is the issue that ETGs with lower stellar 
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Figure 6. Star-formation density (star-formation rate per unit comoving volume) derived by applying the calorimetric technique to the 
far-infrared background radiation (left) and to the optical/near-infrared background (right). In the left-hand panel, the red line shows 
the results for sample 1 and the blue line the results for sample 2. 


masses, such as dwarf elliptical galaxies, do have low values 
of the Sersic index (Graham and Guzman 2003). Allen et al. 
(2015) have used MAGPHYS to estimate the stellar masses 
of the H-GOODS galaxies using the photometry provided in 
Guo et al. (2013). Figure 8 shows the distribution of stellar 
masses for the HGOODS galaxies, which shows that their 
low values of the Sersic index cannot be attributed to low 
stellar masses. 

With a few caveats, the most important of which is 
probably our assumption of a universal initial mass function, 
our conclusion that most stars form in LTGs but most of 
the stellar mass-density today is in ETGs seems secure. A 
necessary implication of our conclusion is that there must 
have been some process that transformed the morphologies 
of many galaxies from the LTG class to the ETG class after 
most of the stars in each galaxy had been formed. 

However, before we consider the mechanisms that might 
cause this major morphological transformation of the galaxy 
population, we have to consider what our LTG/ETG classi¬ 
fication means when we apply it to galaxies at high redshift. 
At low redshift, an LTG is generally a galaxy whose struc¬ 
ture is dominated by a disk and an ETG is a galaxy which 
has a prominent bulge or spheroid. A common assumption 


is that this must also be true at high redshift, which is part 
of the reason for the conclusion from HST imaging stud¬ 
ies that star-forming galaxies at high redshift have struc¬ 
tures dominated by disks (e.g Wuyts et al. 2011; Bruce et 
al. 2012). However, there is at least one way in which galaxy 
evolution might have made this an unsafe assumption. The 
bulges of low-redshift spiral galaxies often have Sersic in¬ 
dices < 2.5 (Balcells et al. 2003). Therefore, one possibility 
we have to consider is that the high-redshift galaxies are 
naked spheroids, which will subsequently accrete a disk and 
grow into the spiral galaxies we see around us today. 

There are two main arguments againsts this idea. The 
first is the circumstantial argument that most of the high- 
redshift galaxies do actually look like they have disks. In 
Figure 9 we have shown the CANDELS images in the Hab- 
band (1.6 (im) of the galaxies in sample 10. There are very 

2 We have chosen sample 1 rather than sample 2, despite the 
fact the latter constitutes statistically a larger fraction of the 
background radiation, because we can be more confident that 
individual galaxies in sample 1 are emitting far-infrared radiation 
since each was detected by the Herschel Space Observatory at 160 
pm. 
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Figure 7. The relationship between the stellar mass-density in the Universe today and the Sersic index of the galaxy in which the 
stars were formed. The two panels show the results for our calorimetric samples of the far-infrared (left) and optical/near-infrared (right) 
backgrounds; the blue and red line in the left-hand panel show the results for samples 1 and 2, respectively, and the blue line in the 
right-hand panel shows the results for sample 3. The histogram for sample 1 in the left-hand panel has been normalised so that the area 
under it is the same as for sample 2. In both panels the black line shows the star-formation rate in the Universe today as a function of 
Sersic index (Fig. 1) normalized so that it encloses the same area as the other histograms. Note the absence of a peak at n = 4, which 
should be present if stars that are currently in elliptical galaxies were born in galaxies with the same structures. 


few galaxies in this sample, even ones with high values of 
the Sersic index, for which there is not some sign of a disk 
on the HST image. 

The second argument that our operationally-defined 
LTGs are similar to the disk-dominated systems in the Uni¬ 
verse today comes from the results of a different technique 
for investigating the morphologies of high-redshift galax¬ 
ies. A number of groups have inferred the morphological 
composition of a high-redshift sample by modelling the dis¬ 
tribution of ellipticities of the images of the galaxies (e.g. 
Ravindranath et al. 2006; Van der Wei et al. 2014). Using 


a large sample of ~40,000 star-forming galaxies from CAN- 
DELS, Van der Wei et al. (2014) show that the distribu¬ 
tion of ellipticities of star-forming CANDELS galaxies with 
1.5 < z < 2.0 and stellar masses > 10 10 Mq can be repro¬ 
duced by a model in which ~75% of the galaxies are disks, 
with almost all the remaining objects having a spheroidal 
structure. Both studies also found a significant population 
with highly elongated structures which could not be ex¬ 
plained by the predicted fraction of highly edge-on disks, 
but in the study of the CANDELS sample these highly elon¬ 
gated galaxies are found below a stellar mass of 10 10 Mq, 
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Figure 8. Estimates of the stellar masses of the galaxies in 
sample 1 from modeling their spectral energy distributions with 
MAGPHYS (Allen et al. in preparation). 


whereas most of the galaxies in our sample 1 have higher 
stellar masses (Fig. 8). 

Despite these arguments, we can not rule out the idea 
that high-redshift LTGs are physically a different type of 
object to the disk-dominated galaxies in the local Universe, 
especially as high-redshift galaxies are physically smaller 
(Bruce et al. 2014) and have more turbulent velocity fields 
(Genzel et al. 2014). However, even if the LTGs at 
high redshift are physically distinct from the disk- 
dominated objects that are LTGs at low redshift, 
this does not eliminate the need to find a mech¬ 
anism for transforming LTGs into ETGs. Suppose, 
for example, the high-redshift LTGs are naked spheroids. 
There still needs to be some mechanism to turn these into 
ETGs, since the obvious evolutionary process, the growth 
of a disk around the spheroid, would still produce a low- 
redshift galaxy with n < 2.5. 

Before we turn to the mechanisms that might have 
caused this morphological shift, it is worth pointing out that 
our result is a statistical result and does not mean there can 
not be evolutionary mechanisms that either leave the mor¬ 
phologies untouched or work in the other direction. For ex¬ 
ample, Graham et al. (2015) have argued that the compact 
passive galaxies discovered at high redshift a decade ago 
(Daddi et al. 2005; Trujillo et al. 2006) are the ancestors of 
the bulges of low-redshift SO or spiral galaxies, and that the 
evolution from one into the other will occur by the growth 
of a disk around the naked bulge - the same mechanism that 
we considered above for the general galaxy population. It is 
possible that this morphological transformation occurs, but 
the following argument implies that its scale is much less 
than the morphological transformation from LTGs to ETGs 
implied by the calorimetry results. 

Graham et al. have estimated that the descendant pop¬ 
ulation in this transformation, the spirals and SOs in the 
Universe today with bulge structures very similar to the 
structures of the compact passive galaxies at high redshift, 
have stellar masses of ~ 10 11 Mg and a volume density of 
3.5 x 1CD 5 Mpc~ 3 dex~ 3 , although this is strictly only a 
lower limit. The stellar mass-density in the descendant pop¬ 


ulation is thus ~ 3.5 x 10 6 Mg Mpc~ 3 . However, the total 
stellar mass-density in the Universe today (the integral of 
the distribution in Fig. 2) is 1.9 x 10 8 M 0 Mpc~ 3 , approxi¬ 
mately 50 times larger. Thus this evolutionary channel may 
exist, but the morphological change from LTGs to ETGs 
revealed by the calorimetry results involves a much larger 
fraction of the current stellar mass-density. 

We now consider the mechanisms that might have pro¬ 
duced this major change, and from now on we assume that 
the LTGs at high redshift are disk-dominated galaxies. Of 
the theoretical ideas for morphological transformation that 
were discussed in §1, two are clearly consistent with the idea 
that LTGs predate ETGs: 1) the idea that ETGs are formed 
by the merger of two LTGs (Toomre 1977); (2) the idea 
that the spheroid in a galaxy is built up by the rapid mo¬ 
tion of star-forming clumps towards the centre of the galaxy 
(Noguchi 1999; Bournaud et al. 2007; Genzel et al. 2011, 
2014). 

The evidence from HST imaging that star-forming 
galaxies at all redshifts are LTGs and passive galaxies are 
ETGs (Bell et al. 2004a; Wuyts et al. 2011; Bell et al. 2012; 
Bruce et al. 2012; Buitrago et al. 2013; Szomoru et al. 2013; 
Tasca et al. 2014) is evidence that the process responsible for 
the morphological transformation is the same one that is re¬ 
sponsible for quenching the star formation in a galaxy (Bell 
et al. 2012). This evidence favours the second of the ideas 
above, since it provides a clear explanation for both morpho¬ 
logical transformation and quenching. Additional circum¬ 
stantial evidence for this conclusion is that most ETGs in 
the Universe today contain some evidence of a residual disk 
once you look closely enough (Krajnovic et al. 2013). 

To explain our result, the morphological transformation 
of an LTG must have occurred after most of the stars in the 
galaxy had formed. The average star-formation rate in the 
Universe dropped rapidly over the redshift range 1 > a > 0 
(Hopkins and Beacom 2006; Eales et al. 2015), suggesting 
that this redshift range may have been the one in which mor¬ 
phological transformation was most active. Two other results 
point in this direction. First, many studies have shown that 
~50% of the stellar mass in ‘passive galaxies’, which over¬ 
laps substantially with ETGs, was formed since z ~ 1 (Chen 
et al. 2003; Bell et al. 2004; Faber et al. 2006). Second, Tasca 
et al. (2014) have observed a gradual decline in the optical 
emission from disks over this redshift range. If this is the 
critical redshift range, the high-resolution wide-held surveys 
with Euclid will be crucial for investigating the transforma¬ 
tion process. 


6 CONCLUSIONS 

In this paper, we have tried to quantify the importance of 
morphological transformation in galaxy evolution. We have 
used the operational definition that an early-type galaxy 
(ETG) is one with a Sersic index > 2.5 and a late-type 
galaxy (LTG) is one with a Sersic index of <2.5. Using this 
definition, we have obtained the following two observational 
results: 

1 We have used the results of the Galaxy and Mass As¬ 
sembly project and the Herschel Astrophysical Terrahertz 
Large-Area Survey to show that in the Universe today 51% 
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of the stellar mass-density is in ETGs but 89% of the rate 
of production of stellar mass-density is occurring in LTGs. 

2 The extragalactic background radiation contains all 
the energy generated by nuclear fusion in stars since the 
Big Bang. We have resolved the background radiation into 
individual galaxies using the deepest far-infrared survey 
with the Herschel Space Observatory and a deep near- 
infrared/optical survey with the Hubble Space Telescope 
(HST). Using measurements of the Sersic index of these 
galaxies derived from the HST images, we estimate that 
~83% of the stellar mass-density formed over the history 
of the Universe occurred in LTGs. 

Our second result is subject to a number of caveats, the 
most importance of which is that our calorimetric method is 
based on the assumption of a universal initial mass function. 
Subject to these caveats, the difference between the fraction 
of the stellar mass-density today that is in LTGs and the 
fraction that was formed in these galaxies implies that there 
must have been a major morphological transformation of 
LTGs into ETGs after the formation of most of the stars. 

Our remaining conclusions are less certain. Since the 
star-formation density started to decline at z — 1 and since 
there is evidence for the build-up of stellar mass-density in 
passive galaxies since that redshift, it seems likely that the 
morphological transformation from LTGs to ETGs occurred 
over the redshift range 0 < 2 < 1. If this conclusion is true, a 
high-resolution investigation of the morphologies of galaxies 
over this redshift range, such as will be possible with Eu¬ 
clid, will be important for determining the process that is 
responsible for the transformation. 

The fact that HST imaging programmes show that star¬ 
forming galaxies at all redshifts are dominated by disks, 
while passive galaxies have spheroidal structures, implies 
there is a single process responsible both for the morpho¬ 
logical transformation and the quenching of galaxies. This 
favours models such as those in which the growth of a bulge 
in the centre of a disk shuts down star formation by stablis- 
ing the disk against gravitational collapse. 
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Figure 9. The H-band images (1.6 fim) of each of the Herschel sources in sample 1. The images are centred on the Herschel positions 
and have a size of 5 x 5 arcsec 2 . The cross shows the near-infrared counterpart to the Herschel source that we have identified using our 
probabalistic analysis. The circle has a radius of 5 arcsec (the remaining sources are shown in the on-line version of the paper). 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 


© 2002 RAS, MNRAS 000, [2] ?? 

















H-ATLAS/GAMA: The Morphological Transformation of Galaxies 25 



Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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Figure 9 - continued 
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